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Abstract 

In  recent  papers  [1,  2]  a  model  for  the  underpotential  deposition  of  Cu  on  Au(lll)  in 
the  presence  of  bisiilfate  ions  was  proposed.  In  this  model  it  was  asstuned  that  the  bisulfate 
ions  formed  a  y/3  x  y/3  template.  This  template  leaves  a  honeycomb  lattice  of  free  sites  for 
the  adsorption  of  copper.  The  clear  implication  is  that  the  first  peak  has  2/3  of  a  monolayer 
of  Cu.  The  second  peak  corresponds  to  the  replacement  of  the  bisulfate  by  copper  in  the 
adlayer.  We  showed  also  that  the  broad  foot  of  the  first  peaJc  is  due  to  a  second  order 
hard  hexagon  like  transition,  which  is  seen  experimentally  by  Itaya[3]  and  Kolb  [4].  The 
interpretation,  based  on  the  STM  and  LEED  observations,  that  the  first  peak  corresponds 
to  only  1/3  of  a  monolayer,  is  consistent  with  our  model  if  it  is  the  bisulfate  ion  that  is 
actually  seen  in  those  experiments.  In  the  present  work  we  discuss  further  refinements  of 
this  theory:  we  show  that  from  the  dimensions  of  the  bisulfate  ion,  an  oxygen  protrudes 
1.84A  above  the  plane  of  copper  adsorption,  emd  the  hydrogen  of  the  bisulfate  protrudes 
2.I3A.  Assuming  the  same  geometry  on  a  clean  Au(lll)  surface  the  bisulfate  layer  should 
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stand  about  4.1  A  above  the  surface.  However,  the  STM  cannot  measure  absolute  heights, 
2md  therefore,  both  cases  appear  as  a  \/3  x  y/3  overlayer.  The  dependence  of  the  peak 
position  with  concentration  is  also  discussed  using  a  recently  developed  theory  of  kinetic 
effects. 


1  Theory 

The  underpotential  deposition  (UPD)  [5,  6]  of  metals  frequently  involves  phase  transi¬ 
tions  that  are  observed  in  the  voltammograms  as  sharp  spikes  [7,  8].  In  earlier  work  [9,  10] 
we  introduced  a  sticky  adsorption  site  model  to  study  phase  transitions  that  occur  when 
the  adsorbate  is  commensurate  with  the  substrate.  This  situation  is  known  to  occur,,  for 
example  in  the  UPD  of  copper  on  gold  (111)  in  the  presence  of  bisulfate  ions.  As  has  been 
shown  in  earlier  work  [1,  2,  11]  the  particular  structure  of  this  voltammogram  is  due  to  a 
sequence  of  first  and  second  order  transitions  taking  place  on  the  surface.  To  obtain  sharp 
spikes  in  a  voltammogram  the  substrate  surface  must  be  a  perfect  single  crystal  with  large 
domain  size.  Slow  scanning  rates  are  also  required.  This  immediately  suggests  that  thermo¬ 
dynamic  equilibrium  and  reversibility  are  necessary  for  the  occurence  of  these  spikes,  which 
are  associated  with  first  order  phase  transitions  in  the  adsorbate  layer  [9,  10]. 

We  should  mention  in  this  context  the  early  work  of  Bewick  and  Thomas  [12]  who  studied 
the  case  of  UPD  of  Pb  on  Ag  (111),  and  the  work  of  Buess-Hermann[13]  on  the  adsorption 
of  alcohols  on  mercury.  In  both  cases  the  mean  field  Frumkin’s  isotherm  [14]  was  used  to 
interpret  the  observations.  We  observe  here  that  in  both  cases  the  adsorbed  layer  is  incom¬ 
mensurate  with  the  substrate,  and  in  our  theory  the  adsorbate  must  be  commensurate.  For 
the  commensurate  case  Frumkin’s  isotherm  corresponds  to  the  mean  field  theory,  which  is 
known  to  be  very  inaccurate.  The  coexistence  curve  is  of  the  form  y  =  for  Frumkin, 
while  the  exact  lattice  result  is  y  =  x^/*  [15].  This  has  important  consequences  for  the  shape 
of  the  curves  and  for  the  sensitivity  of  the  voltammogram  to  the  interaction  parameter  of  the 
adsorbate.  In  the  case  of  the  exact  lattice  gas  isotherm  the  occupancy  9cu  changes  suddenly 
from  an  almost  Langmuir  type  behaviour  to  a  phase  transition  at  the  critical  value  of  the 
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interaction  parameter  g-i-  In  the  mean  field  theory  the  change  is  much  more  gradual. 

For  the  incommensurate  case,  or  for  the  liquid  mercury  electrode,  Frumldn’s  isotherm 
describes  a  gas-liquid  transition,  and  not  the  solid-liquid  transition  [16].  This  opens  inter¬ 
esting  possibilities  about  different  types  of  trainsitions  that  are  possible,  and  that  we  hope 
to  discuss  in  the  future.  In  the  voltammogram  the  area  under  the  spike  is  proportional  to 
the  charge  transferred,  not  to  the  coverage  of  the  surface,  because  the  charge  per  adatom  on 
the  surface  is  not  necessarily  an  entire  number  equal  to  the  stoichiometric  electrovalence  [17]. 

In  this  paper  we  discuss  further  refinements  and  consequences  of  our  model  [1,  2].  We 
analyze  the  specific  geometry  predicted  by  it,  we  present  a  detailed  discussion  of  the  depen¬ 
dence  of  the  model  parameters  on  the  potential  bias,  and  we  discuss  briefly  am  extension 
to  the  djmamic  case  for  slow  voltage  scanning  rates.  We  use  these  results  to  explain  re^nt 
experiments  in  which  the  concentration  dependence  was  studied. 


2  THE  MODEL 

The  case  of  the  UPD  of  Cu  on  Au(lll)  in  the  presence  of  H2SOa  has  been  extensively 
studied  in  recent  times  both  experimentally  [3,  18,  19,  20,  21,  22,  23,  24,  25,  26,  27]  and 
theoretically  [28].  We  should  also  mention  that  two  recent  studies  of  the  kinetics  of  this 
system  [29,  30]  clearly  indicate  that  in  the  shape  of  the  voltammogram  the  first  peak  can  be 
derived  from  equilibrium  considerations,  as  was  done  in  our  earlier  work,  while  in  the  second 
the  kinetics  is  much  slower  and  therefore  observable  from  the  ^malysis  of  the  voltammogram. 


In  recent  papers  [1,  2]  a  model  for  the  underpotential  deposition  of  Cu  on  Au(lll) 
in  the  presence  of  bisulfate  ions  was  proposed.  We  summarize  here  the  main  results  of  our 
previous  work. 

In  our  model  for  the  underpotential  deposition  of  Cu  on  Au(lll)  in  the  presence  of 
sulfuric  acid,  we  assume  that  a  well  defined  sequence  of  events  takes  place: 
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•  Bisulfate  ions  are  adsorbed  at  very  positive  potential,  forming  a  \/3  x  \/3  lattice  on  the 
gold  surface.  These  bisulfate  ions  are  desorbed  as  the  potential  is  decreased,  undergoing 
a  hard  hexagon  like  second  order  phase  transition. 

•  Copper  ions  are  then  adsorbed  on  the  free  adsorption  sites.  The  adsorption  of  copper 
produces  a  reabsorption  of  the  bisuUate  ions,  which  eventually  will  undergo  the  hard 
hexagon  transition  to  rebuild  the  \/3  x  lattice  on  the  gold  surface.  This  forms  a 
honeycomb  lattice  for  the  adsorption  of  the  remainder  of  the  copper  of  the  first  peak. 

•  In  the  second  peak  the  adsorbed  bisnlfate  ions  are  displaced  by  copper  ions  from  the 
y/Z  X  >/3  positions.  However,  they  could  still  remain  bound  to  the  copper,  which  now 
forms  a  full  monolayer  on  top  of  the  Au(lll)  surface. 

We  assume  that  the  bisulfate  binds  to  the  gold  (111)  surface  in  such  a  manner  that  the  sulfur 
is  directly  on  top  of  the  adsorption  site  for  the  copper,  three  of  the  bisulfate  oxygens  b^g 
above  and  directly  associated  with  the  three  gold  atoms  of  the  surface,  which  form  a  triangle 
about  the  adsorption  site.  Packing  considerations  indicate  that  two  bisnlfate  groups  caxmot 
be  adsorbed  onto  neighboring  adsorption  sites.  The  bisulfate  ions  will  thus  form  a  y/3  x 
film  by  occupying  one  of  the  three  triangular  sublattices  At  of  the  full  triangular  lattice  of 
adsorption  sites,  with  a  maximum  coverage  of  1/3. 


3  BISULFATE  ADSORPTION 

In  our  model  it  was  assumed  that  the  bisulfate  ion  formed  a  \/3  x  y/d  template.  This 
template  leaves  a  honeycomb  lattice  of  free  sites  for  the  adsorption  of  copper.  The  clear  im¬ 
plication  is  that  the  first  peak  has  2/3  of  a  monolayer  of  Cu.  The  second  peak  corresponds 
to  the  replacement  of  the  bisulfate  by  copper  in  the  adlayer.  We  showed  also  that  the  broad 
foot  of  the  first  peak  is  due  to  a  second  order  hard  hexagon  like  transition,  which  is  seen 
experimentally  by  Itaya[3]  and  Kolb  [4].  We  believe  that  the  interpretation  that  the  first 
peak  correponds  to  only  1/3  of  a  monolayer,  based  on  the  STM  and  LEED  observations,  is 
consistent  with  our  model  if  it  is  the  bisulfate  ion  that  is  seen. 
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In  the  present  conununication  we  amplify  the  previously  proposed  HB  model  [1,  2].  In 
this  theory  the  broad  foot  of  the  first  spike  in  the  Cu- Au  voltammogram  is  due  a  second  order 
surface  pheise  transition,  similar  to  the  hard  hexagon  phase  transition  [31],  This  transition 
was  actually  observed  by  Itaya  [3]  and  by  Kolb  [4]. 

In  our  theory  in  the  initial  stages  of  the  process  there  is  a  strong  coadsorption  of  copper 
with  the  bisulfate.  At  positive  potentials  {V  >  0.45  volts  with  respect  to  the  standard 
(Ag/AgCl)  electrode),  the  bisulfate  is  strongly  adsorbed  onto  the  clean  Au(lll)  surface. 
This  means  that  there  is  an  ordered  \/3  x  \/3  structure  for  these  very  positive  potentiab 
observable  by  STM  or  AFM.  We  assume  that  it  retains  its  charge,  and  therefore,  the  bisulfate- 
bisulfate  interaction  is  both  long  ranged  and  repulsive.  The  HSO^  sits  in  a  tripod  position, 
that  is  with  its  three  oxygen  atoms  almost  directly  atop  the  Au  atoms  of  the  substrate,  so 
that  the  adsorption  of  one  HSO^  necessarily  excludes  nearest  neighbor  occupation.  Xbis 
makes  the  short  ranged  part  of  the  surface  interaction  mathematically  isomorphic  to  the 
hard  hexagon  problem,  solved  mathematically  by  Baxter[31].  There  will  be  a  second  order, 
order-disorder  phase  transition  when 


es>ec  =  0.2764,  (1) 

where  6s  is  the  fraction  of  the  Au(lll)  lattice  adsorption  sites  that  are  occupied  by  the 
bisulfate. 

Consider  a  model  in  which  there  are  three  species,  which  for  brevity  we  will  call  E  (empty 
sites),  Cu( copper)  and  S  (  bisulfate).  In  our  lattice  model  of  the  surface  there  is  no  inter¬ 
action  between  E  and  the  other  adsorbates,  S  strongly  repels  S,  and  S-Cu  and  Cu-Cu  are 
strongly  attractive.  Models  in  which  three  components  are  adsorbed  were  recently  discussed 
in  the  literature  [32,  33,  34,  35]. 

We  define  [9,  10]  the  inner  layer  equivalent  fugacity  z$  for  the  adsorption  of  the  bisulfate 

zs  =  XsWPsiO,rk)  (2) 

where  the  sticking  coefficient  can  be  interpreted  as  XsH’)  =  exp[/3^5],  with  fis  as  the  free 
energy  change  that  occurs  when  a  bisulfate  ion  binds  to  the  gold  surface,  removing  in  the 
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process  the  adsorbed  water.  is  the  inner  layer  local  density  of  bisulfate  for  a  local 

potential  0. 

We  assume  the  simple  exponential  form  [36]  for  the  fngadty 


where  0  =  If  kT  is  the  Boltzmann  thermal  factor,  the  electrosorption  valency  of  the  bisulfate 
is  ^5  =  —1)  a^d  rl;§^  is  the  electrosorption  reference  potential,  that  depends  on  the  nature 
of  the  substrate. 

The  adsorption  isotherm  for  the  hard  hexagon  model  has  been  derived  by  Joyce  [37].  His 
results  can  be  fitted  to  Fade  approximants  [2].  For  the  high  density,  ordered  phase 

\2 


,,,,  0.2764  +  0.155(z5  -  11-09)  +  0.01(^5  -  11-09)^ 

^  1  +  0.5(25  -  11.09)  +  0.03(zs  -  11.09)2 


(4) 


For  the  low  density,  disordered  phase, 


_  0.70925  +  0.00792|  +  O.OOllzl 


0lOW  ^ 


0  <  25  <  11.09. 


(5) 


1  4-  3.325 

Consider  now  the  adsorption  of  bisulfate  onto  Au(lll)  in  the  absence  of  copper.  Using 
Eq.(3)  with  the  effective  charge  of  the  adsorbed  bisulfate  ^5  =  —1, 


^5  = 


(6) 


where  we  take 


tpp  =  0.388^  (7) 

and  the  bisulfate  fugadty  25  is  (2^  =  1,  T  =  298.16iif) 

Z5  =  e38-922tt^-^#“]^  (8) 

The  bisulfate  is  desorbed  as  the  potential  turns  less  positive.  With  these  parameters  the 
S  lattice  becomes  disordered  at  about  0.45  V,  which  is  when  the  copper  starts  to  deposit. 
As  long  as  the  S  template  is  disordered,  the  Cu  is  adsorbed  randomly  (probably  next  to  the 
bisulfate). 

The  copper  will  start  depositing  at  around  ^  =  0.4V,  but  since  the  bisulfate  is  in  its  dis¬ 
ordered  phase,  there  should  be  no  cooperative  effects  due  to  the  copper-copper  attraction  on 
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the  surface,  and  therefore,  to  a  first  approximation  we  may  assume  that  the  electrodeposition 
of  copper  follows  a  Langmuir-like  (  or  Frumldn-Langmuir  [36])  adsorption  isotherm 


ecu  =  Bs  <  =  0.2764  (9) 

where  V’cu  **  the  electrodeposition  reference  potential  of  Cu  on  gold  (111),  and  Ccu  is  the 
effective  electrovalence  (certainly  related  to  the  electrosorption  valency  )  of  Cu  (see  the 
discussion  below,  Eq.(27)). 

The  electrosorption  of  bisulfate  on  polycrystalline  copper  [38,  39]  shows  that  it  binds 
much  more  strongly  than  to  gold.  We  may  assume  that  the  electrosorption  potential  of  the 
bisulfate  varies  linearly  with  the  amount  of  copper  in  the  surface.  Thus,  guided  by  the  work 
of  Trasatti  on  the  influence  of  the  anion  on  the  potential  of  zero  charge  [40],  we  take 


(10) 

In  our  model,  then 

_  -0  -Csi9e(t^-t/’s) 

*5  —  ■*56  , 

(11) 

with 

i>S  =  +  (1  -  ^Cu)V’^", 

(12) 

In  the  presence  of  copper  the  bisulfate  is  then  readsorbed.  Eventually,  the  fraction  of 
occupied  sites  will  again  surpass  the  critical  value  of  0.2764  (which  corresponds  to  83%  of  the 
occupancy  of  the  template  >/3  x  y/Z  lattice),  and  then  there  will  be  a  first  order  transition 
for  the  copper  on  the  honeycomb  lattice  of  the  Cu  sites. [41]  If  we  assume  that  the  occupancy 
of  copper  Bcu  =  1/3  at  the  transition,  we  get  that  i/’f”  =  -.302V,  which  is  in  qualitative 
agreement  with  Wieckowski.  [38,  39]  In  figure  1  we  show  the  adsorption  isotherm  of  the 
bisulfate  as  a  function  of  potential  [2]. 


4  THE  FUGACITY  OF  COPPER 

We  recall  that  in  the  using  the  sticky  site  model  [9, 10],  there  are  two  basic  parameters 
that  determine  the  phase  behaviour  of  the  adsorbed  layer: 
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•  The  loc2d  fugadty  z  =  zcu  of  the  two  dimensional  adsorbate  gas  is  the  product  of  the 
sticking  probability  Acu  and  the  contact  density  Pcu(0»V’)  of  the  ion  t, 

ic«  =  Ac„pL(0,V’)  (13) 

•  The  lateral  interaction  parameter  ^2 

^2=52(Ri,Rj)  [i^j  =  nearest  lattice  neighbors],  (14) 

where  R^,  Rj  are  the  positions  of  the  neighboring  lattice  sites  i,j,  and  g^ir)  is  the  pair 
correlation  function  in  the  electrode  plane,  but  in  the  absence  of  the  adsorption  sites. 

Both  of  these  parameters  depend  on  the  appUed  potential  rj^.  Consider  first  the  fugadty: 
The  sticking  probabihty  is  a  function  the  bonding  free  energy  ecu  of  the  adsorbate  to  the 
binding  site.  This  quantity  is  independent  of  the  ionic  strength  of  the  electrolyte,  for  a  gjyen 
potential.  However,  it  will  depend  on  the  potential  bias  [36]  at  the  surface.  We  write 

Ac.  =  Aj.e'’*'^-,  (15) 

where  the  binding  energy  of  the  Cu  ion  depends  on  the  potential  bias  of  the  metal  substrate, 
or  electrode, 

ecu  =  CcuefV’  -  ^^u]y  (16) 

where  e  is  the  elementary  charge,  V*  is  the  potential  at  the  electrode  surface,  and  is  the 
electrosorption  reference  potential,  (cu  is  the  effective  charge  of  the  adsorbed  spedes  at  the 
surface,  not  necessarily  equal  to  the  electrovalence  in  solution. 

The  contact  density  Pc:ij(0)  V’)  is  a  also  function  of  the  potential  bias.  As  a  first  approxi¬ 
mation  we  wiU  aissume  that  the  contact  density  is  given  by 

/cuM  =  (17) 

where  /3c;u(0,0)  is  the  contact  density  at  zero  potential  xl>,  which  for  the  Gouy-Chapman 
theory  is  equal  to  the  bulk  density  pcuy  and  evcu  is  the  full  charge  of  the  Cu  ion.  Therefore 

^cu  =  ^cuPcuiO^  (18) 
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Similarly,  we  write  for  the  lateral  interaction  parameter  gi 


92  =  (19) 

where  uj2  is  the  potential  of  mean  force  between  moeties  adsorbed  in  neighboring  sites.  Since 
we  know  that  the  charges  of  these  moeties  change  with  V’[42],  ijJ2  must  also  be  a  function 
of  the  potential  bias  (  A  detailed  theory  of  a  model  which  calcTilates  this  effect  will  be 
presented  separately.) 


Therefore  we  may  expand  about 

w,  =  +  ... 


(20) 


where 


a 


1 


du)2 


As  was  discussed  in  earlier  work  [9,  10],  the  coexistence  curve  is  given  by  the  exact 
condition 

^  =  1,  (22) 

where  qi  is  the  number  of  neighbors  of  the  lattice,  6  for  the  triangular  (center  filled  regular 
hexagons)  lattice  and  3  for  the  honeycomb  (  empty  center  hexagons).  If  the  parameter  u 
is  less  than  unity  then  we  are  in  the  1  phase  region.  If  it  is  larger  than  1,  then  we  are  in 
the  two  phase  region.  The  coexistence  curve  is  obtained  setting  u  =  1  in  Eq.(22).  For  the 
triangular  lattice  [43] 


Ot  =  (1/2) 


16flf2 


1)  (^2  +  3). 


with  a  critical  value  =  Z.  For  the  honeycomb  lattice  [41] 


Oh  =  (1/2) 


V  (95"-i)’(SJ-1)7  ^ 


with  a  critical  value  =  (2  +  \/3)^  —  13.93. 


(23) 


(24) 
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As  it  was  shown  in  earlier  work,  [10]  the  adsorption  isotherm  can  be  expressed  as  a 
generalized  Langmuir  adsorption  isotherm  of  the  form 

a  _  >^(^1  “) 

i  +  A{z,uy 

where  u)  =  ^  is  a  polynomial  in  z  and  u.  For  the  case  in  which  there  is  a  phase 
transition,  9cu  is  practically  a  step  function  which  was  represented  by  an  error  function  in 
our  previous  work.  A  very  simple  and  useful  alternative  form  is 

A  =  z  +  {g2-l)u-,  (26) 

where  n  is  some  entire  number.  The  width  of  the  transition  spike  is  inversely  proportional 
to  n  (actually  to  n/4)  ,  and  will  play  the  same  role  as  the  width  parameter  of  our  previous 
work[44].  It  is  related  to  instrumental  width,  substrate  domain  size  and  defects  on  the  sur¬ 
face.  - 

Combining  Eq.(18)-  Eq.(20)  and  Eq,(22),  we  get  that  at  the  transition  point,V’T 

The  relation  of  the  parameters  in  this  formula  to  the  electron  transfer  coefficients  will  be 
discussed  in  future  work. 


5  KINETIC  MODEL 

There  is  a  very  extensive  literature  on  the  voltammograms  for  the  electrodeposition  of 
metals,  starting  with  the  early  work  of  Koryta  [45],  Delahay  [46]  Brdicka  [47],  Guidelli[48], 
Laviron  [49],  Lorenz[50],  Reinmuth[51],  Sadkowski[52]  and  Lantelme[54].  However  these 
models  do  not  correspond  to  the  case  of  adsorption  onto  a  lattice  with  cooperative  effects. 
In  our  general  statistical  mechanics  treatment  of  phase  transitions  in  the  ad-layers  [10]  the 
adsorption  isotherm  is  written  in  the  form  of  a  Fade  approximant,  which  represents  to  any 
required  precision  the  exact  coexistence  line  and  universal  critical  behaviour.  This  form  of 
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the  isotherm  lends  itself  to  a  simple  way  to  include  kinetic  effects,  which  in  the  limit  of  very 
slow  scanning  rates  gives  back  the  required  equilibrium  phase  behaviour. 


A  kinetic  equation  for  the  fraction  of  occupied  sites  6cu  of  the  honeycomb  or  triangular 
sublattices  of  the  original  lattice  can  be  obtained  for  a  simple  probabilistic  model,  in  which 
the  probability  of  a  site  to  be  occupied  by  Cu  is  proportional  to  the  density  of  Cu  at  the 
electrode  surface  and  the  number  of  free  sites. 


/>Ca(0,  V')(l  —  Bcu) 

The  probability  of  an  ion  that  is  adsorbed  to  be  dissolved  is  proportional  to  9cu-  Therefore, 
the  rate  of  change  is 

^  =  KaPcuiO.  0)(1  -Bcu)-  KaBcu  (28) 

The  two  reaction  rate  constants  Ka  and  Kj  must  satisfy  at  equilibrium  the  relation 
Eq.(25)  and  therefore 

KapCuM)  ^  (29) 

Kd 

There  are  many  ways  in  which  this  relation  can  be  satisfied:  If  Kd  is  of  the  exponential  form 


1361 

Kd  = 

(30) 

then 

i^a/»c„(O,0)  =  ^foe“^M(u), 

(31) 

We  write  Eq.(28)  in 

the  form 

^  =  -ecuK,d  +  K,pcu{0,rj^), 

(32) 

where 

Kad  =  KaPCu{0,^)  +  Kd. 

(33) 

This  is  our  nonlinear  kinetic  equation,  which  together  vdth  the  diffusion  equation  consti¬ 
tutes  a  system  of  equations  that  can  be  solved  numerically.  Interesting  asymptotic  solutions 
for  limiting  cases  can  be  obtained,  that  urill  be  discussed  in  a  separate  paper. 


If  we  neglect  diffusion  then  an  asymptotic  solution  for  slow  scanning  rates  can  be  obtained. 
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We  consider  here  the  case  in  which  0  is  the  electrode  potential  and 

0  =  0m  +  vi  (34) 

where  0in  is  the  initial  value  of  the  potential,  and  v  is  the  scanning  rate. 

The  solution  of  Eq(32)  is  [55]: 


Bcuit)  =  +  [9cu{0)  -  <?^„(0)]  e-/o  (35) 


where 


/loo  /A\  _  ^aii)PCu{Q',t)  _  •4.(ucu(f))  /o-\ 

~  l  +  A{ucu{t)y 

is  the  equilibrium  fraction  of  occupied  sites  for  the  potential  at  time  t.  Changing  now 
variables  to  0,  and  integrating  by  parts  we  get  the  expansion  in  terms  of  v  :  i  - 


[Kaiirk)  50 


(37) 


KadW  50  ■  KM)d^ 

This  equation  shoTild  be  useful  in  calculating  the  overall  kinetic  constants  Kad  bom  the  shifts 
in  the  voltammogram. 


A  simple  useful  formula  that  can  be  obtained  from  this  equation  is  that  for  small  scanning 
rates  t;  the  zero  scanning  rate  limit  is  given  by 

«?.W  =  (l/2)[«Sf'W  +  «ciW].  (38) 

where  B^^{tp)  is  the  isotherm  for  oxidation  scan  and  Bq^{iI^)  for  the  reduction  scan. 


6  THE  CURRENT  IN  THE  VOLTAMMOGRAM 

Consider  the  electrode  interface,  in  which  a  current  flows.  We  have  specifically  in  mind 
the  case  of  the  UPD  of  Cu  onto  Au(lll).  In  this  case  there  are  two  ions  that  participate  in 
the  transport  of  charge: 
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•  The  bisiilfate,  that  carries  a  negative  charge  which  it  keeps  when  adsorbed  onto  the 
surface.  Therefore,  it  will  contribute  only  to  the  capacitive  current  density  jc-  The 
other  contribution  to  jc  is  from  the  diffuse  layer. 

•  The  copper,  which  is  adsorbed  and  discharged  to  some  extent.  It’s  electrovalence  will 
change  from  ucu  in  the  bulk  phase  to  (cu  at  the  electrode  interface.  Furthermore,  (cu 
should  be  a  function  of  the  potential  i/>. 

Therefore,  the  total  current  density  is  (  for  a  recent  discussion  see,  for  example  the  work  of 
De  Levie  [53]  and  Lantelme  [54]) 

JT  =  jc  +  jp  (39) 

where  jj  is  the  total  current  density,  jc  is  the  capacitive  contribution  amd  jp  is  the  Faradaic 
contribution.  If  we  neglect  double  layer  effects,  then  the  capacitive  current  is  due  to^he 
bisulfate,  and  the  faradaic  current  almost  exclusively  to  the  discharge  of  the  copper  ions. 
Then 

ic  =  (lM)[G  +  ^f]|  (40) 

gives  the  current  due  to  the  discharge  of  the  capacitor  as  well  as  its  change  in  integral  ca¬ 
pacitance  Ci.  A  is  the  area  of  the  electrode. 

The  current  associated  with  the  cation  (in  this  case  the  copper  )  can  be  written  as 

j>  =  (A//A)e  -ec.~  +  (vo.-Cc.)^  ^  (41) 

where  M  is  the  number  of  adsorption  sites  per  area  A,  e  is  the  elementary  charge,  (cu  is  the 
actual  charge  of  the  adsorbed  copper  amd  t{>  the  potential. 

In  accordance  with  the  quantum  theory  [42],  the  charge  of  the  adsorbate  will  change  with 
the  applied  potential  because  the  electron  density  at  the  surface  is  changing.  Intuitively  we 
take  the  exponential  form 

=  (42) 

where  and  V’l*'  aje  adjustable  parameters. 


13 


In  figure  2  we  show  a  comparison  of  the  current  density  obtained  fi;om  the  copper  con¬ 
tribution  of  our  model  to  the  oxidation  branch  of  Kolb’s  voltammogram.  We  have  adjusted 
the  position  of  the  hard  hexagon  transition  to  .22V,  the  position  of  the  first  peak  to  .21V. 


7  GEOMETRY  OF  THE  MODEL 

We  discuss  now  some  geometrical  properties  of  our  model.  We  recall  that  we  assume 
that  the  bisulfate  ion  formed  a  y/3  x  template.  This  template  leaves  a  honeycomb  lattice 
of  firee  sites  for  the  adsorption  of  copper.  The  clear  implication  is  that  the  first  peak  has  2/3 
of  a  monolayer  of  Cu.  In  this  section  we  want  to  show  that  the  bisulfate  actually  protrudes 
from  the  surface,  and  therefore  the  1/3  of  a  monolayer  seen  by  the  STM  and  LEED  obser¬ 
vations  is  consistent  with  our  model.  ;  . 

The  vertical  distance  between  the  substrate  Au(lll)  plane  and  the  adsorbed  copper  can 
be  estimated  to  be  2.22  A  ,  by  taking  the  mean  distance  in  metallic  gold  and  in  metallic 
copper.  From  early  EXAFS  experiments  [18, 19,  20]  and  recent  theory  [28]  we  know  that  the 
gold  (or  copper  )  distance  to  the  oxygens  in  the  bisulfate  is  between  1.95A  and  2.1  A  ,  so 
that  from  the  dimensions  of  the  bisulfate  ions,  it  should  be  between  1.84A  and,  if  we  count 
the  hydrogen  of  the  bisulfate,  2. 13 A  higher  than  the  plane  where  the  copper  lies.  Assum¬ 
ing  the  same  geometry  on  a  clean  Au(lll)  surface  the  bisulfate  layer  should  stand  about 
4.06-4.35  A  above  the  Au  surface.  However,  the  STM  cannot  measure  absolute  heights, 
and  therefore,  both  cases  appear  as  a  y/Z  x  y/Z  overlayer. 

An  observation  to  be  made  is  that  the  Cu  —  0  distances  appear  at  1.81  A  if  we  assume 
an  undistorted  bisulfate  ion  and  copper  honeycomb  lattice.  This  is  a  bit  lower  than  the 
expected  value  of  1.95 A  .  Therefore  we  would  expect  that  the  adsorbed  ion  be  deformed  so 
that  the  oxygens  lie  about  .2  A  below  the  line  of  the  copper  atoms,  which  would  relieve 
the  stress,  but  at  the  same  time  create  a  activation  energy  barrier  for  the  desorption  of  the 
copper.  This  again  is  consistent  with  the  experimental  observations. 

The  proposed  coordinates  are  shown  in  Table  1,  and  the  geometry  in  figure  3. 
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8  CONCENTRATION  DEPENDENCE 


If  all  the  parameters  stay  the  same  as  the  concentration  of  copper  is  changed,  a  shift 
in  the  concentration  will  necessarily  imply  a  shift  in  the  transition  potential 

kT 


-Aln(/)cu(0,0)]. 


e{vcu  +  9La/2  -  Ccu) 

For  a  tenfold  increase  in  concentration  the  potential  shift  is 

59.2mV 


A^t  = 


where  !/«//  is  the  effective  charge  obtained  directly  from  the  equation  Eq.(27) 

I'eff  =  t'Cu  +  qLOt/2  —  Cctt- 


(43) 


(44) 


(45) 


We  are  now  in  position  to  discui^s  the  recent  experiments  of  Omar  et  al.  [30]  and  Holtzle 
et  al  [29].  Using  the  results  of  the  local  kinetic  theory  [55]  we  can  estimate  the  positions 
of  the  peaks  in  the  experiment  of  Omar  et  al.  We  notice  that  the  relative  shift  of  the  bulk 
deposition  edge  corresponds  exactly  to  the  electrovalence  of  copper,  2.  We  should  remark 
that  this  a  bit  surprising,  since  one  would  expect  that  it  should  scale  to  the  activity  of  the 
copper,  not  its  concentration.  However  this  also  means  that  the  contact  density  pcu(0,V’) 
in  our  model  scales  almost  exactly  to  the  Gouy-  Chapman  estimate  Eq.(17),  which  is  non¬ 
trivial  and  reassuring.  We  recall  that  our  interface  is  not  planar,  our  solution  is  molecular, 
and  the  position  at  which  pcu(0,  V')  is  taken  is  some  average  over  the  positions  in  the  metal 
electrolyte  interface.  We  know  that  the  contact  theorems  are  satisfied  by  these  averages,  and 
this  is  probably  the  reason  why  the  Gouy-Chapman  estimate  is  so  accurate. 


If  we  use  the  shifts  of  the  edge  to  compute  the  effective  charge  Eq.(45)  for  the  peaks  in 
the  voltammogram  we  find  that  the  effective  charge  for  the  second  peak  is  20  %  higher  than 
2,  and  that  for  the  first  peak  we  get  about  1,  which  is  very  difficult  to  explain  in  terms  of 
individual  fugadties  or  the  Nemst  equation  [30]. 
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The  inclusion  of  the  variation  of  the  pair  correlation  functions  between  adsorbed  particles 
with  ^  removes  these  difficulties.  The  results  of  table  2  show  that  the  remarkable  accuracy 
of  the  Gouy-Chapman  estimate  holds  also  for  other  transition  points.  However  the  actual 
values  of  the  charges  cannot  be  reconciled  mth  the  other  experimental  evidence. 

For  the  second  peak  a  charge  bigger  than  two  cannot  be  explained  in  any  other  way  but 
by  including  the  variation  of  with  potential  The  value  q  =  0.14  seems  re2isouable. 

The  apparent  charge  for  the  first  peak  Veff  =  1  is  definitely  inconsistent  with  the  exper¬ 
imental  evidence  and  in  particular  with  the  integrated  charge  shown  in  fig.6  by  the  same 
authors,  where  the  ratio  of  the  charges  of  both  peaks  is  very  closely  1,  so  that  each  peak  cor¬ 
responds  to  50%  of  the  monolayer  charge.  If  we  assume  only  1/3  coverage  by  copper  ais  has 
been  indicated  in  the  literature  [3, 4],  the  charge  would  be  only  16%  of  the  monolayer,  cleuly 
impossible.  The  explanation  that  the  remainder  of  the  charge  is  due  to  desorption  of  the 
sulfate  is  inconsistent  with  the  EXAFS  experiment  [18],  with  the  radiotracer  experiment  [38], 
and  with  the  microbalance  experiment [26],  all  of  which  indicate  that  there  is  strong  bisulfate 
adsorption  even  after  the  full  monolayer  is  deposited. 

The  situation  improves  if  we  take  the  2/3  monolayer  coverage  predicted  by  the  HB 
model  [1,  2).  We  get  50  %  of  the  monolayer  for  each  peak’s  charge  when  i^cu  —  Ccu  —  10 
and  Ccu  =  10  .  The  value  ^Cu  =  0.5,  which  implies  i/g//  —  Ccu  =  1-5  satisfies  the  60%  for 
the  first  peak  charge  requirement  of  Kolb’s  votammogramm,  and  produces  a  neutral  film 
near  the  expected  value  of  the  point  of  zero  charge  for  this  system.  We  believe  that  this  last 
one  is  a  more  plausible  value  of  the  charge,  since  charge  neutrality  explains  the  lifting  of  the 
recunstruction  of  the  gold  surface  by  the  adsorbed  copper  layer. 
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Table  1 


STRUCTURE  OF  THE  OVERLAYER 


Caption  Table  1 

Structure  of  the  adlayer  with  1/3  coverage  of  bisulfate  and  2/3  of  Cu.  The  stressed  structure 
was  made  using  an  undistorted  bisulfate  ion.  In  the  relaxed  structure  the  bisulfate  ion  was 
squeezed  so  that  the  0  —  0  distance  is  shortened  from  2.42  A  to  2.30  A  ,  while  the  copper- 
oxygen  goes  from  1.81  to  1.87  A. 
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Table  2 

SHIFTS  OF  THE  VOLTAMMOGRAM’S  PEAKS 


Exp. (mV) 

Shifts(mV) 

e 

#1 

#2 

#3 

#2-#l 

#3-#2 

Eff.Chg. 

Real.Chg. 

a 

Edge 

299 

270 

241 

29 

29 

2.04 

0 

Foot 

557 

500 

458 

57 

42(?) 

- 

1st.  peaJc 

554 

496 

444 

58 

52 

1.02/1.14 

0.5 

-0.33/-0.25 

2st.  peak 

352 

327 

303 

25 

24 

2.37/2.47 

0 

0.14 

Caption  Table  2 

Shifts  of  the  voltammogram’s  peaks  with  concentration  from  Omar  et  al  [30],  Figure  4. 
Experiment  #  1:  9.10~^ff2SO4,  5.10“*C'u5O4 

Experiment  #  2:  ^.1Q~^CuSOa  -  ■ 

Experiment  #  3  :  9.10~^/f2-5O4,  5.10“‘‘Cu5O4 

The  effective  charge  is  evaluated  using  Eq.(45).  a  is  evaluated  using  Cc«  =  0.5  for  the  first 
peak,  and  (cu  =  0  for  the  second  peak. 

The  foot  corresponds  to  the  hard  hexagon  transition,  which  we  cannot  estimate  accurately 
from  the  graphs. 
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FIGURE  CAPTIONS 


FIGURE  1 

Fraction  of  sites  occupied  by  HSO^.  Dotted  line  represents  the  clean  gold  surface.  Solid 
line  the  case  with  copper  coadsorption.  This  graph  is  qualitative  [2]. 


FIGURE  2 

Comparison  of  the  theory  to  the  oxidation  branch  of  Kolb’s  voltammogram  [21],  for 
O.liV  H2SO4,  1Q~3N  CuSO^  and  a  scamning  rate  of  ImVfs.  The  theory  includes 
only  the  copper  contribution. 


FIGURE  3 

Geometry  of  the  adlayer  with  1/3  bisulfate  and  2/3  of  copper.  The  copper  ring  holds  the 
bisuUate  structure,  but  the  top  of  the  bisuUate  is  still  higher  than  the  copper. 
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